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Abstract 

We examine the zero-range limit of the finite square well in arbitrary dimensions through a systematic 
analysis of the reduced, s-wave two-body time-independent Schrodinger equation. A natural consequence 
of our investigation is the requirement of a delta-function multiplied by a regularization operator to model 
the zero-range limit of the finite-square well when the dimensionality is greater than one. The case of two 
dimensions turns out to be surprisingly subtle, and needs to be treated separately from all other dimensions. 
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I. INTRODUCTION 



Every undergraduate student of physics familiarizes themselves with the mathematical ma- 
chinery of quantum mechanics by solving the one-body, one-dimensional (ID) time-independent 
Schrodinger equation (TISE) for a variety of potentials which admit an exact solution. One need 
only glance at any number of introductory quantum mechanics textbooks 1 to note that the attrac- 
tive finite square well (FSW) and delta function potentials are among such standard potentials. 
These potentials provide a pedagogical introduction to the concepts of bound and scattering states 
in ID, while illustrating some of the pathologies that can occur when the potential is singular at a 
point. 

Owing to the nature of the boundary conditions imposed on the wave function by the one- 
body potential, the mathematical analysis required for the solution of the bound and scattering 
states of the ID FSW and delta-potential are quite different. However, many textbooks make 
a point of commenting that the ID delta potential results can be recovered from the ID FSW 
of depth Vo and range b, in the limits V — >■ oo, and b — > 0, provided that the area under the 
well, V b, remain constant (see e.g., problem 2.31 in Ref. [ l|]). A question then naturally arises: 
can higher dimensional zero-range interactions be similarly constructed by taking an analogous 
limiting procedure? 

Zero-range interactions have received renewed attention in the context of harmonically trapped 
cold atoms,- where the dilute nature of the gas and ultra-low temperatures permit a treatment of 
the interactions between the charge neutral atoms as being primarily two-body and of zero-range, 
in the relative s-state. The symmetry of the relative s-state allows for a reduction of the two- 
body problem to an effective one-body problem that is already familiar to undergraduate students 
exposed to elementary scattering theory. Given that the dimensionality of these systems can be 
manipulated in the laboratory, 3 theoretical studies of such interacting, ultra-cold atomic systems 
involve the use of zero-range interactions appropriately generalized to two and three dimensions 
(2D and 3D, respectively). 

In this paper, we present a systematic extension of the analysis in ID for obtaining a zero-range 
interaction from a FSW to arbitrary dimensions. Our goal is to illustrate, using only elemen- 
tary quantum mechanics, how zero-range interactions in higher dimensions are constructed by 
starting from a well-studied finite range potential. We will focus solely on two-body interactions 
in the relative s-state (I = 0), with a view to readers interested in applications to, e.g., modern 
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ultra-cold atoms research. We point out that the naive replacement of a finite square well in the 
6—7-0 limit with a bare delta-potential is only valid in ID, whereas in higher dimensions, the bare 
delta-potential must be "regularized" to avoid mathematical divergencies. Our regularization of 
the delta-potential provides an intuitive introduction to the notion of pseudo-potentials in the con- 
text of zero-range interactions without having to invoke the more technical language of Green's 
functions and self-adjoint extensions that are sometimes used in the literature. 

The organization of our paper is as follows. In the next section, a brief introduction the hyper- 
spherical coordinates is presented, which then sets the stage for our general analysis of the d- 
dimensional FSW in Sec. III. In Sec. IV, we examine the b — > limit of the FSW in arbitrary 
dimensions. In Sec. V, we present our concluding remarks. 



H. HYPER-SPHERICAL COORDINATES 

In the following analysis, we make use of rf-dimensional hyper-spherical coordinates. 4 These 
coordinates are comprised of a radius, r, and d — 1 angular coordinates. In 3D, for example, the 
coordinates would be r and the two angular variables 9 and <p. In what follows, we consider only 
the relative coordinate {i.e., the centre-of-mass, R = (r + r')/2, has already been removed so that 
we have an effectively one-body problem), and as such, the <i-dimensional hyper-radius is defined 
as 



a 



\ i=i 



where Xi is the i-th component of the first particles position, r. Owing to our focus on the relative 
s-state two-body interaction, the two-body wave function of the system will be hyper-spherically 
symmetric. It is useful to note that in this case, the integration of a function f(r) over a hyper- 
sphere of radius R is given by 

dir d / 2 r R 

/CO*" = FTTTT^T / fir)r d - x dr , (2) 



IW2 + 1, Jo 

where the pre-factor is the rf-dimensional analogue of the familiar 3D "47r" result. The use of the 
hyper-radius is straightforward in all dimensions except perhaps ID, where some confusion may 
arise. Following the definition in Eq. CO the ID hyper-radius is 



\[{x — x') 2 = \x\. (3) 
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Thus, while r G (0, oo), the Cartesian variable x G (—00, 00). Finally, note that in ID the pre- 
factor of the integral given by Eq. (2) is 2, which reflects the fact that the analogous ID s-state 
function is even so that 

/R pR pR 

f(x)dx = 2 f(x)dx = 2 f(r)dr . (4) 
-R, JO Jo 

III. THE FINITE SQUARE WELL IN ARBITRARY DIMENSIONS 

We now consider the d-dimensional generalization of the FSW, which takes the form 

V(r) = -V Q(b - r) (5) 

where Vq > is the strength of the interaction and b is the range. The rf-dimensional TISE for 
/ = then reads 

-^V><(r)-^<(r) = ^<(r) r <b (6) 

and 

h 2 

-^V^>(r) = ^>(r) r>b (7) 

where M is the mass of each particle, ^<(r) and ^>(r) denote interior and exterior solutions, 
respectively, and = + 



A. Scattering States: E > 



1. Interior solution: r < b 



Let us first consider the region r < b (the so-called interaction region). With rf = M(Vq + 

E) /h 2 = MV /h 2 + k 2 we can write Eq. © as 

c^<(r) d-ldjj^r) 2 

, 2 + j = -r]^<{r). (8) 

dr z r dr 

This ordinary differential equation (ODE) is solved^ by a linear combination of Bessel functions 
of the first, J a , and second, Y a , kind with a = d/2 — 1. However, since the r < b region includes 
the origin, we discard Y a as a possible solution as it diverges at the origin. We then write our wave 
function in the interior region as 

Cl 7]b d l 2 J a {r]r) 

where c\ is a constant of integration. 



2. Exterior solution: r > b 



Next, we consider the asymptotically free exterior region (r > b), where the TISE is given by 

+ d-lijp^ = 
dr z r dr 

with k 2 = ME/h 2 . The solution to the above ODE can be written as 5 

giJ a (kr) + g 2 Y a (kr) 

^>( r ) = ^72=1 > ( n ) 

where g± and #2 are constants of integration. Note that here, we retain the Y a solution, in contrast 
to the interior region where it had to be ignored. 



3. Boundary conditions at r = b 

Given that our FSW potential has no singular behaviour, we require that the logarithmic deriva 
tive of the wave function be continuous at the boundary r = b 



^<( r )Jr=b \^>( r )/r=b 



(12) 



where ' = d/dr. 

The result of Eq. (12) is 



, 9iJ a +i{kb) + g 2 Y a+1 (kb) \ = J a +i(rjb) 

giUkb) + g 2 Y a (kb) J V U V b) ' ^ ) 



which allows us to write the ratio 21 as 
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£fr _ kJ a+1 (kb)J a (r]b) - r]J a (kb)J a+1 (r]b) 
gi r]J a+1 (r]b)Y a (kb) - kY a+1 (kb)J a (rjb) ' 

We may then re- write Eq. (1 1) in the form familar from 3D scattering theory, 7 viz., 



(14) 



Mr) = A iUkr) ~X-i° Yaikr)) ' (15) 

where A is a constant which is generally dimensionally dependent, and 5 is called the s-wave 
phase shift.™ The s-wave phase shift characterizes the strength of the scattering in the I = partial 
wave by the potential V(r), at the energy E = h 2 k 2 /M. Explicitly, we have 

(x \ - V J a(kb)J a+ i(T]b) - kJ a+1 {kb)J a {rjb) 
anl 0j " V J a+1 ( V b)Y a (kb) - kY a+1 (kb)J a ( V b) • 1 } 



Now, as k — >■ 0, the short-distance physics should become irrelevant, and the scattering should 
not depend on the details of the short-range two-body potential. With this in mind, we expand 
Eq. (16) in the k — > limit, retaining only the leading order term in k,- viz., 



7T I k\ 1 



— tan($o) — — ttt~\ — m r ~ 7 r, d ^ 2 ( 17 ) 

yo) T(d/2)T(d/2-l)\2j b 2-d 1 + 2^ Mf) \ 

where now, 77 = V / Mp ! .^ We observe that for d ^ 2, we may define the following quantity, 

a d -* = ^ 1 lim tan((5o) = 1 fl8) 

7T r(tZ/2)r(d/2 - 1) fc-vo A; d " 2 ft2-a ^ I 2-d JaW x 

^ T)b J a+1 (r)b) ^ 

from which we may re- write Eq. (17) as 



7T 



d-2 



ta °W- r W 2)r W 2-i) iYj ' <19) 

The crucial point here is that the s-wave phase shift, Eq. (19), now depends on only a single 
quantity, a, and is independent of the shape of the potential. The quantity a then completely 
characterizes the k — > limit of the scattering, so that any short-range potential having the same 
value of a will lead to absolutely identical scattering. The single parameter a is called the s-wave 
scattering length, which is now generalized to d ^ 2 by Equation (18). It is easy to see that if we 
put d = 3 in Eq. (17), and apply Eq. (18), we obtain an explicit connection between the scattering 
length and the range, b, of the FSW 

, = ,(,.!!«), (20) 

which is well-known for the 3D finite square well. For completeness, it's also worthwhile noting 
that for d = 1, we obtain 

^/MV /h 2 b 



which agrees with the symmetric ID result discussed in great detail in Ref. 
The expression analogous to Eq. (17) for d = 2 is given by 



J2U 



tan5 ~ — 7777^' rf = 2 ' 

rjbJx (rjb) 



2 (lnf + 7 + ln6 + 

where 7 ~ .577215665... is the Euler constant. In this case, the 2D s-wave scattering length is 
defined by 



2g2tan(5 T J Q (r,b) 

lim = be^-h(vb) (23) 

fc->o k 



and Eq. (22) becomes 

- tan<5 ^ - — - * a — r , (24) 
2 (In f + 7) 

again illustrating that in the k — > limit, the scattering is completely characterized by the s-wave 
scattering length, a. 

The rather formal definition of the scattering length, viz., Eqs. (18) and (23), for d 7^ 2 and 
d = 2, respsectively, can be given a more familiar geometric interpretation as follows. Let us first 
recall the low energy behaviour of J a and Y a £ 

J ° (kr) - nm (t) • <25) 



and 



, - (ln(^) +7) d = 2 

7r \kr) ' ' 



from which we obtain the A; — s- scattering solutions, viz., 

tan<5 r(d/2- l)r(d/2) / 2 



V>>(r)~|l+ u w ^ (^) h rf ^ 2 ' ( 27 ) 



d-2 N 



^>(r)~ 1 ° In— +7 , d = 2. (28) 



and 

. . / 2t 
ib^ it I ~ I 1 

77 V 2 

Inserting Eq. (19) into Eq. (27) and Eq. (24) into Eq. (28), we obtain (to within an overall unim- 
portant constant) 

^>(r)~l-^, c^2, (29) 

and 

^>--^> 

Let us first discuss a > (attractive potential). In this case, the s-wave scattering length is 
the value at which the k — > scattering wave function obtains a node (for r > 6). Thus, the 
low-energy scattering wave function is "bent over" so that it has a negative slope at the boundary, 
implying that it can be joined to an exponentially decaying solution corresponding to a genuine 
bound state. Therefore, the k — > scattering state can be connected with a shallow bound state 
when a > 0. Conversely, a < (repulsive potential), indicates the absence of a bound state. 
Indeed, it is the backward extrapolation of the k — > wave function to r < which provides 
a measure of the negative scattering length. An excellent, and more thorough, discussion of the 
interpretation of the s-wave scattering length and its sign can be found in Ref . Ill ill . 



It is interesting to note that in 2D, the scattering length, Eq. (23), is always positive. Within 
this context, note that in ID and 2D, a > no matter how small the depth, Vo, of the well. This 
observation leads to the "well-known" conclusion that in ID and 2D, an arbitrarily weak attractive 
short range interaction will always support at least one bound state. 12 Loosely speaking, in ID and 
2D, any depth, V , is enough to bend the wave function so that a > 0, whereas in 3D, a minimum 
depth is required before the first bound state occurs; for the 3D FSW, the minimum depth is 
readily seen to be Vq > 7i 2 h 2 /4Mb 2 .— In the very special case of 2D, where the scattering length 
is always positive (or a — > oo when Ji in Eq. (23) has a zero), there are only bound and loosely 
bound states^- for an attractive short-range potential. Indeed, a — > oo signals the appearance of the 
next loosely bound state, which develops into a bound state as the depth of the potential is further 
increased. In the next subsection, we will develop the formal connection between scattering states, 
with a > 0, and bound states. 

B. Bound States: -V Q < E < 

We have already briefly introduced the notion of bound states supported by the FSW in any 
dimension. In this section, we will focus our attention on energies for which — Vq < E < 0, 
and thereby make the connection between low-energy bound states and the low-energy scattering 
states discussed above. In order to accomplish this goal, we must address how the condition 
— Vo < E < affects our solutions to the TISE in the asymptotically free (r > b) region; the 
solutions will be unaffected in the interior region since rf > still remains true. 

1. Exterior solution: r > b 

With E negative, k = ME / H 2 is imaginary, and we can write k = in where k = 
y^—ME/h 2 . We may then use our scattering solution, viz., Eq. (11), with k — > in. The iden- 
tification k — > %k is only a formal tool, and if desired, a direct solution of Eq. (10) with — k 2 
replaced by k 2 could be pursued instead. 

To begin, we make use of the following properties of Bessel functions of imaginary argument: 5 



J a {inr) 



(31) 



and 



Y a (iKr) 



i a+1 I a (Kr) - -i~ a K a (Kr) , 



(32) 
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where I a and K a are the modified Bessel functions of the first and second kind, respectively. 
Utilizing these expressions in Eq. (11) allows us to write for the exterior solution, 

iMr) = ' 3/a ^i 4/aM , (33) 

where g% and are dimensionally dependent integration constants. Of the two modified Bessel 
functions, I a increases exponentially, and as a result we set g 4 = 0. Thus our exterior, normalizable 
bound state wave function reads^^ 

*>«■) = (34) 

2. Interior solution: r < b 

As mentioned earlier, the interior solution is again given by Eq. (9) owing to the fact that for 
-Vo < E < 0, rj 2 > still holds true. 

3. Boundary conditions at r = b 

The allowed energies of these bound states are determined by matching the wave function and 
its derivative at the boundary r = b. As a result, the energy is determined by solutions to the 
equation 

><(r)/ r=6 \^>( r )Jr=b 

which is identical to Eq. (12), but now applied to bound states. In ID, Eq. (35) reproduces the 
standard result for the even bound states 1 

k — Tj tanrjb. (36) 

For d-dimensions, Eq. (35) reads 

r K a+1 (Kb) = ^ J a+ i(yb) ^ 
K a (nb) J a {r]b) 

Our intention is to now relate the low energy bound states to the discussion of the scattering length 
given in the previous subsection. To proceed, we put k = %k in Eq. (16), use Eqs. (31) and (32), 
along with Eq. (37), to obtain the remarkably simple result in all dimensions: 

cot 5q = i . (38) 
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Therefore, the bound states can be directly obtained from the E > scattering states in any di- 
mension provided we put k — >• in and cot(5o) = i- Recall that for scattering states, the continuity 
of the logarithmic derivative of the wave function at r = b fixed the s-wave phase shift. In the 
present case, the same boundary condition at r = b has likewise fixed the phase shift, but now, 
the phase shift is a purely imaginary number cot(5 ) = i, associated with the purely imaginary 
momentum, k = in. We wish to point out that in standard treatments of scattering theory, the same 
result, viz., cot(5 ) = i, is obtained, but involves the introduction of the 5 (fc) -matrix (or equiva- 
lently, the partial-wave scattering amplitude), and an analysis of its analytic properties.-^ Here, 
we have accomplished the same goal, but without having to introduce any additional mathematical 
machinery. While our approach may not be as mathematically elegant, it requires less formalism, 
and is therefore more accessible to students with only a limited exposure to scattering theory. 

Using Eq. (38), we may finally make the connection between the low-energy bound state en- 
ergy, E = —h 2 K 2 /M and the s-wave scattering length. Quite simply, Eqs. (19) and (24) are 
evaluated at cot(<5 ) = i, giving 



and 



2e~ 7 

k = i , d = 2. (40) 

a 



Since k = Imag[/c], we obtain in both ID and 3D the well known result 1 *^ k = 1/a, and E 
—h 2 /Ma 2 for the shallow bound state energies. In 2D, we obtain k = 2e _7 /a giving E 
—Ah 2 e~ 21 /Ma 2 for the shallow bound state energy.— 



IV. CONTACT INTERACTION LIMIT OF THE FINITE SQUARE WELL 

Following the standard ID treatment of allowing a FSW to go to a zero-range interaction (i.e., 
delta-function) we wish to investigate the b — > limit of the above developed results. For simplic- 
ity, our analysis will be formulated in terms of the low-energy bound states, although exactly the 
same results will also hold true for the low-energy scattering states; this is not surprising in view 
of our discussions up to now. We will take the b — > limit of the FSW while insisting that the 
(i-dimensional area inside the well remain constant. This area is given by 

V = / V{r) dr = - rr , /9 , U V . (41) 

All Space r(d/2 + 1) 
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The negative sign in Eq. (BP) merely reflects the fact that our <i-dimensional FSW is attractive and 
thus V(r) < while r > 0. To keep Vq constant, we require that 

illustrating that as & — >■ 0, V — > oo. We may now write our TISE in terms of V . We obtain 

fr* 2l/ , V T(d/2 + l)e(b-r) ,, , , . 

where it is to be understood that, in place of ijj{r) above, we use ij} < for r < 6 and ^> for r > b. A 
standard analysis of a one-dimensional attractive delta function interaction centered at the origin^ 
involves integrating the TISE in a neighbourhood about the origin from — e to e, followed by 
letting e — > 0. We extend this approach to the case of arbitrary dimensions in the sense that we 
now integrate our above TISE over a hyper-sphere of radius e > b, and then allow b — > (this 
gives us a contact interaction) followed by e — > 0. We will focus on the d ^ 2 case now, and 
present the 2D result at the end of this section. 

We recall that for the bound states, we have the following solutions: 

Cl T]b d l 2 J a (rjr) 

Mr) = < 44 > 

and 

Mr) - ^ ■ (45) 
The Laplacian of the exterior solution is given by 



2 cib d / 2 rf J d /2+i(yr) _ dbWrfdJdpjrir) 

' J d /2(vb) r"/ 2 - 1 J d/2 ( V b) r<*/ 2 



Integrating the TISE over a hyper-sphere of radius e > b we obtain 

(rr=b pr=e 
/ ^<(r)dr+ / ^>(r)dr 
</r=0 Jr=fe 

When we let 6 — > 0, followed by e — > 0, it is straightforward to see that in this limiting procedure, 

-r=b rr=b 



f V^<(r)dr, f V<(0^, rV>(r)dr, (48) 

Jr=0 Jr=0 7r=6 



will all vanish. 
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s 9j rw^i) i + (49) 



Next, anticipating 6 — >■ 0, we expand ip>(r) for small r and obtain (for d ^2) 

1/2 - 1) /2\ d/2_1 J_ 

An application of the d-dimensional generalization of the divergence theorem, v/z., 4 ' 16 ' 17 

\r d -y r(d/2 + i) 1 J ' 1 ; 

gives us 
The integral 

w/2+i) r^<M, T=Ci ^ o (52) 



7r d / 2 ./ r=0 & d 



is independent of 6. We may then write Eq. (1471) as 

^ 3 r(rf/2-l)c/(d-2)7r d / 2 /2\ d/2_1 



2Mr(d/2 + l) W 

from which we find 



cidVb = , (53) 



S3 VT{d/2-l) fJ n , d/2 f2\ d/2 - 1 



V ^'i 2MT(d/2 + ^ d - 2 ^\-.) ■ (54) 

This form for V ensures that our TISE has the proper behaviour under integration in the 6—^0 
limit. This expression for V gives us a potential that reads 



ci 2M v V K / b^o ¥ 
The issue now is that the potential in Eq. (55) only acts on ip < (r), but in the zero-range limit, 
ip < (r) will be "squeezed out" as the entire interior region is reduced to a point. Mathematically, 
the remedy to this situation is to insist that 

^ 9(6-r) = f AaS^MO^^Mdr , (56) 

where O^- 1 is an operator to be determined, and is a dimensionally dependent constant. is 
sometimes referred to as a regularization operator ,— and its neccessity can be traced back to the 
fact that \/d ^ 1, ip> (r) is undefined as r — > 0. Thus, if we were to try to use a bare delta function 
in the integral on the right-hand side of Eq. (56), the integral would be ill-defined \fd^l. 
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The integral on the left-hand side of Eq. (56) is 



hm — ^— ; — '-ibJridr = — =-. , (57) 

b^o b d ' r (d 2 + 1 ' v ' 



whereas for the right-hand side we have 



f Aa*M(r)0M^>(r)rfr = A^Jm6« (^^) • (58) 
Equating the results of Eqs. ([57]) and (l58l) . we obtain the constant A d , viz., 

A Cl 1 

We are now free to replace lim^ 6> ^d" r ' > with Ad5^(r)0^ in Eq. (55) to get, for the 6 — > limit 
of the FSW, 

= _ft^M (J _ (2\«" *<4(r) , M (60) 



2Mr(d/2 + i) 'W Um „ 6«(^e 

The appropriate operator, 0^ d \ has different forms depending on the dimensionality of the system. 
Specifically, is found by requiring that the singular behaviour of lim r ^ ^>( r ) m Eq. (56) is 
removed. Let us orient ourselves first with ID, and then move on to higher dimensional spaces. 
In ID, ^>(r) is regular as r — y 0, so the proper operator is = 1 and we have 

= ft 2 ^ g%) (61) 

M lim^o ( v/f AL 1/2 («r) ) 



We note that lim r ^ (V^^~i/2( Kr )) = -\/ 7r/2^, and we obtain the known result 

V(r) = -^W(r) = -^«(r), (62) 

where we have made use of k = 1/a for the shallow bound state found in Sec. Ill B above. 
Equation (62) is to be viewed as the zero-range interaction reproducing the same A; — ?• scattering 
as for any short-range potential with the same scattering length, a. Again, note that in ID, the 
6—7-0 limit of the FSW is proportional to a bare delta function potential. 

For higher dimensions (d ^ 2), the operator needed is = ( 2 r d ~ 2 , with 



?™\dr) r^ 2 " 1 2r(d/2) \2J ' K } 
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In defining the operator 0^ d \ we have kept to the common convention in the literature and used the 
partial derivative, d/dr,to emphasize that OW only acts on the radial component of any function 
it encounters. Of course, if acts on a function of r only, the partial derivative is to be treated 
as a full derivative. 

The expression for the b — > limit of the FSW then becomes 

(-i)^r(d/2)r(d/2- 1) (2 



_ h 2 dn d / 2 
^ ~ Afr(d/2 + l)r(d-2) 

For example, setting d = 3 in Eq. (64) gives 



7T \ K 



d-2 



d\ d — 2 



*«(r)(£J 

(64) 



= ^~ 5 (r) ¥ r ' (65) 

where again, for low-energy bound states, k = I /a. Equation (65) is in perfect agreement with 
what is found using other approaches in the literature.-^^ - — 

For 2D, we need to step a little bit further back. The small r expression given in Eq. (|49| ) is 
valid for all d ^ 2. For <i = 2 we have 

^> (r) ^ -9s (In nr/2 + 7) , (66) 

with 

V^>(r) = -27r^5 (2) (r) . (67) 

Except for the term involving Vl^^r), a ^ °f me integrals in Eq. (|47|) are the same in 2D as they 

are in any other dimension. As a result, our integrated TISE for 2D is 

n2g 3 h 2 
M 

and we obtain 



+ 2 Cl V = , (68) 



So in 2D we have 



V(r) = - 9 ^rn^±. (70) 

Again, we wish to replace the part involving the Heaviside function with an operator involving the 
delta function. By the same argument as for all other dimensions, we obtain 

A 2 = ^ ^ , (71) 

93 linw)0( 2 ) (#o(«r)) 
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and our 2D zero-range potential becomes 



V(r) = ¥® o<, (?2) 

Af lim r _„0<2> (#„(«•)) 



The proper operator for 2D is 0^ 2 ^ = 1 — r In 2e ^_ 7 J^, and we note that 



lim I 1 -rln-^— 4" ) #o(«0 = -1 . (73) 



r^o y 2e 1 "T <9r 

It follows that the zero-range limit of the FSW in 2D is given by 

V(r) = S ( ' (r) 1 - r In — ; — — 

2vr/l 2 



M 



5(2) ( r ) ( l-rln— 4" I > ( 74 ) 
\ aeor J 



where Eq. (40) has been used to connect the low-energy bound state to a scattering state with 
a > 0. 



V. CONCLUSIONS 



We have presented a systematic analysis of the FSW in arbitrary dimensions, thereby providing 
generalizations to quantities typically introduced in the context of 3D scattering; namely, the s- 
wave scattering length and phase shift . 7 We have shown that the d-dimensional scattering length, 
a, can be naturally interpreted as the node (or its extrapolation) of the k — > wave function. In the 
k — > limit, we have also illustrated how the s-wave scattering length completely characterizes 
the scattering properties of the system for any short-range potential in any dimension. For students 
familiar familiar with 3D scattering, the 2D results of the FSW highlight that even in simple "toy- 
model potentials", dimensionality plays a pivotal role in determining the physical properties of the 
system. For example, we point out that in ID and 2D, a > for any Vo, which introduces the 
notion that an arbitrarily weak short-range attractive interaction in ID and 2D will always support 
at least one bound state .— 

In order to avoid introducing additional mathematical formalism, we have also shown how to 
utilize only the scattering states with imaginary momentum, in, to obtain imaginary phase shifts, 
which for a bound state near E — > 0~, can be connected to the s-wave scattering length discussed 
earlier. This alternative approach has the virtue of utilizing mathematical techniques familiar from 
the study of elementary ID problems i.e., continuity of the wave function and its deriative at the 
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boundary to treat both bound and scattering states in a unified way. We feel that students may 
benefit from this presentation, particularly those not yet exposed to the analytic properties of the 
S(k) -matrix jililS 

Following our treatment of the d-dimensional FSW, we examined its b — > limit by extending 
the ID analysis^- to arbitrary dimensions. Our main result is that the zero-range limit of the FSW 
is given by 

™ 2MIW2 + 1) { " 2) U UmMO o M (^ ' <?5) 



for d 7^ 2, and 



7( r ) = ^!^(r)fl-rln-^|-) , (76) 
K J M v 7 V 2e 1 -T <9r ' 



for d = 2, respectively. The non-trivial expressions for the b — > limit of the FSW serve to 
illustrate that the usual textbook suggestion of obtaining the delta function results from the limit 
of a FSW is not so straightforward in dimensions greater than one. In fact, the d = 2 case 
turns out to be the most interesting, with subtle mathematical issues not typically discussed in 
undergraduate quantum mechanics courses. Nevertheless, we have tried to use only the TISE and 
a minimal amount of mathematical machinery to motivate the so-called regularized delta function 
potentials given by Equations (75) and (76). In more technical papers, similar results are developed 
in the context of self-adjoint extensions, and Green's function techniques and are given the name 
"pseudo-potentials" . 13 ' 15 ! 19 — 

Finally, it is our hope that the presentation used in this paper may serve as a basis for introducing 
a more general treatment of the FSW in undergraduate quantum mechanics. To this end, we 
suggest the following useful exercise. Initially, the student would be asked to explore Eq. (37) 
using standard graphical solutions, for d — 1 and d = 3, with fixed b, and varying Vo. Next 
the student would be asked to fix Vob d and let Vo — > oo and b — > and also graphically look 
for solutions. The student will find that in the latter limit, the ID version of Eq. (37) has only 
one, finite energy, bound state, while the 3D solution obtains more and more bound states as Vo 
increases. The point of this exercise would be to show that a naive replacement of the ID FSW 
with a delta function potential leads to only one bound state in ID, which is perfectly consistent 
with known results. 23 However, the analogous V — > oo, b — > limit in 3D (with V b 3 fixed) leads 
to an infinite number of bound states, with an unphysical, infinitely bound ground state energy. 
These findings would then motivate a different approach to the zero-range limit of the FSW in 
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3D, namely, the approach suggested in this paper. The student could then be guided to obtain the 
appropriately regularized 3D delta potential, and if desired, the connection between the regularized 
3D delta potential and self-adjoint Hamiltonians^^ could be pursued. 
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